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ABSTRACT 

The Sloan Digital Sky Survey (SDSS) and photometric/spectroscopic surveys in the GOODS-South 
field (the Chandra Deep Field-South, CDF-S) are used to construct volume-limited, stellar mass- 
selected samples of galaxies at redshifts < 2 < 1. The CDF-S sample at 0.6 < z < 1.0 contains 207 
galaxies complete down to M = 4x 10^° Mq (for a "diet" Salpeter initial mass function), corresponding 
to a luminosity limit for red galaxies of M-q — —20.1. The SDSS sample at 0.020 < z < 0.045 contains 
2003 galaxies down to the same mass limit, which corresponds to Mb = —19.3 for red galaxies. 
Morphologies are determined with an automated method, using the Sersic parameter n and a measure 
of the residual from the model fits, called "bumpiness," to distinguish different morphologies. These 
classifications are verified with visual classifications. In agreement with previous studies, 65% — 70% 
of the galaxies are located on the red sequence, both at z ~ 0.03 and at z ~ 0.8. Similarly, 65% — 70% 
of the galaxies have n > 2.5. The fraction of E-^SO galaxies is 43% ± 3% at z - 0.03 and 48% ± 7% 
at z ^ 0.8; i.e., it has not changed significantly since z ~ 0.8. When combined with recent results for 
cluster galaxies in the same redshift range, we find that the morphology-density relation for galaxies 
more massive than 0.5Af* has remained constant since at least z ~ 0.8. This implies that galaxies 
evolve in mass, morphology, and density such that the morphology-density relation does not change. 
In particular, the decline of star formation activity and the accompanying increase in the stellar mass 
density of red galaxies since z ~ 1 must happen without large changes in the early-type galaxy fraction 
in a given environment. 

Subject headings: galaxies: clusters: general — galaxies: elliptical and lenticular, cD — galaxies: 
evolution — galaxies: formation — galaxies: fundamental parameters — galaxies: 
general — galaxies: photometry 



INTRODUCTION 
The morphology-density relation (MDR; iDressled 



Il980l ). i.e., the observation that high-density environ- 
ments contain a higher fraction of early-type galaxies 
than low-density environments, provides a clue about the 
formation and evolution of galaxies, and suggests that en- 
vironment plays an important role in shaping the galaxy 
population. The interpretation of the MDR, however, 
is not straightforward because of additional correlations 
between morphology, galaxy mass, color, star formation 
history, and metallicity; all quantities tha t depend on en- 
vironment (e.g.,.H oe g et al . 2003: Kauffma nn et al.ll2003l . 
120041: iBlanton et al.[l2005at iBaldrv et al.ll2006l ). There 
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are strong indications that the stellar mass of a galaxy, 
independently of th e environment, det ermines the color 
(iBaldrv et al.l2006[ ) and concentration (| Kauffmann et ahl 
|2004[ ) of a galaxv. In addition, the environment plays a 
role in star formation activity (|Kauffmann et aLll2004D . 

It is a matter of debate whether the fate of a galaxy is 
ultimately determined by its mass or its environment. 
Studies of galaxy properties such as those mentioned 
above at higher redshifts help to answer those ques- 
tions, and there is abundant evidence for both scenarios. 
The massive galaxy populations of z ~ 1 clusters ap- 
pear to be fully assern bled and passively evolving (e.g., 
iDe Propris et al.ll2007[ ). whereas the galaxy population 
as a whole undergoes significant evolution between z = 1 
and the present, both in terms of mass and color (e.g.. 
Bell et al] l2004bD an d in terms of star formation (e.g. 



Le Floc'h et al.ll2005l ). This points to a strong relation 



between formation epoch and environment. On the other 
hand, high-mass galaxies have had lower specific star for- 
mati on rates than low-mass galaxies since at least z ^ 2 
(e.g.. lJuneau et all [20051 iNoeske et al.ll2007t l. and there 
is only a very modest age difference between field and 
cluster early-type galaxies with m asses M > 10^^ Mq at 
z ~ 1 (van Dokku m fc van der M arcl 2007), with possi- 
bly a stronger age dependence on ma,ss tha n environment 
(|van der Wei et al.ll2005l: iTreu et al.ll2"005[ ). All these re- 
sults, at first sight contradictory, will have to be recon- 
ciled with each other. 

The MDR, and its evolution with redshift, plays a cru- 
cial role in our understanding of the effect of the environ- 
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ment on galaxy evolution, and significant evolution of the 
early-type galaxy fract ion has been found between z ^ 1 
and t he present fe.g.. iDressler et all 119971 : ISmith et al.1 
[200a iPostman et al.il2005D . These studies are based on 
samples that are selected by luminosity, and to ensure 
that the z ~ 1 samples contain the progenitors of the 
galaxies in the local samples, the luminosity limit for 
samples at different redshifts is corrected for evolution. 
It is assumed that the evolution is the same for all galax- 
ies, which might not be the case. In addition, luminos- 
ity is very sensitive to bursts of star formation that are 
likely far more prevalent at z 1 than in the local uni- 
verse. Therefore, it is worthwhile to adopt an alterna- 
tive selection method that is less sensitive to potentially 
rapidly changing galaxy properties such as the star for- 
mation rate. Stellar mass is a quantity that changes less 
rapidly than luminosity, such that it might be a more 
suitable tracer of the z ^ 1 progenitors of local galax- 
ies. In addition, it is more straightforward to compare 
stellar masses to model predictions. Obviously, mass is 
not an ideal tracer of the evolving galaxy population, 
just like luminosity, as it increases through s tar forma- 
tion and mergers. StiU, in lHolden et"aI1 (|2007l . , hereafter 
Paper I) we show that there is no discernible evolution 
between z — 0.83 and the present in the morphological 
mix of galaxies that are more massive than 4 x 10^° Mq 
(~25% of the mass of a typical cluster galaxy, ~50% of 
the mass of a typical field galaxy). The evolution in the 
early-type fraction seen in luminosity-selected samples 
of cluster galaxies is due to the relative increase in the 
number of blue, star-forming, low-mass galaxies. 

Given the abundant evidence for significant evolution 
since z ~ 1 in the galaxy population outside clusters, it 
is important to measure the early-type galaxy fraction 
at lower densities. It is not inconceivable that the early- 
type galaxy fraction docs evolve strongly in those envi- 
ronments because of the strong evolution of the galaxy 
population in the field and in groups. On the other hand, 
if morphological evolution is primarily driven by envi- 
ronmental effects, the morphological mix of the galaxy 
population in a given density might not evolve. Several 
authors have stud ied the morphological mix of galax- 
ies at z ~ 1 (e.g., 
iCaoak et al. l mm . 



iBeh et al.l[2004alMBTmdv"eraIII^^ 
Abraham et~al] l2007| ). A range of 



techniques and sample selection methods is used in those 
investigations, and it is not always clear how the high-z 
results compare to local studies and how sample selection 
affects the interpretation. In particular, mass estimates 
and morphologies are generally not obtained with the 
same method for local and distant galaxy samples, such 
that systematic errors may contribute to the observed 
evolution in, for example, the early-type fraction. 

In this paper we construct mass-selected samples of 
low- and high-redshift galaxies with morphologies clas- 
sified in an internally consistent manner, such that the 
morphological mix of the galaxy population at low and 
high redshift can be directly compared both with each 
other and with our recent results on samples of cluster 
galaxies over the same redshift range (Paper I). In § 2 
we describe our galaxy samples extracted from the Sloan 
Digital Sky Survey (SDSS) at redshifts 0.020 < z < 0.045 
and the Chandra Deep Field-South (CDF-S) at redshifts 
0.6 < z < 1.0 and determine galaxy masses and mor- 
phologies. In § 3 we present and discuss our results 



and their implications for the evolution of the early- 
type galaxy fraction in our field samples and, through 
comparison with cluster galaxies, how the MDR evolves 
with redshift. Finally, in § 4 we put our results in 
the context of the evolving red sequence and the de- 
cline in the average cosmic star formation rate den- 
sity. We adopt the standard cosmological parameters, 
(r!M, r!A, h) = (0.3, 0.7, 0.7). 

2. DATA 

First, we construct a complete, stellar-mass-selected, 
volume-limited sample of galaxies at redshifts 0.02 < 
z < 0.045 froni the Sloan Digital Sky Survey 
(SDSS: lYork et all l2000l). D ata Release 5 (DR5; 
lAdel man-McCarthv et al.lr2007| ). Similarly, we construct 
a stellar mass-selected, volume-limited sample at 0.6 < 
z < 1.0 in the Chandra Deep Field-South (CDF-S). Then 
we determine the morphologies of the galaxies in both 
samples in a consistent manner and quantify the environ- 
ment (parameterized by the local surface galaxy density) 
in which the galaxies are situated. 

2.1. A Mass-Selected, Volume- Limited SDSS Galaxy 

Sample 

The sample consists of two redshift bins: galaxies at 
redshifts 0.02 < z < 0.03 with stellar masses 4 x 10^" < 
M/Mq < 1.6 X 10", and galaxies at redshifts 0.035 < 
z < 0.045 with stellar masses M > 1.25 x 10" Mq. 
The reason for this split is to ensure spectroscopic com- 
pleteness on both the faint and the bright end and to 
minimize systematic problems with the photometry of 
bright galaxies (see below). The g-band magnitudes of 
the galaxies in these samples range from gmod = 13.7 
to 16.5. Here (/mod is the model magnitude from the 
SDSS pipeline (see below for the definition). The up- 
per limit (7inod = 13.7 ensures that this sample is not 
incomplete due to the cutoff at bright magnitude s in the 
SDSS spectroscopic survey (jStrauss et al.l[200^ . This 
source of incompleteness is essentially a surface bright- 
ness limit imposed to avoid saturation. The lower mag- 
nitude limit (7mod = 16.5 ensures spectrosc opic complete- 
ness at the faint end (jBlanton et al.|[2005bD and the feasi- 
bility of determining morphologies without being limited 
by background noise. Finally, only those galaxies located 
at least 130 pixels from the edge of the image tiles are 
included in the sample, in order to avoid problems with 
the morphological analysis described below. This does 
not introduce systematic effects other than decreasing 
the volume of the sample by 28%. The lower mass limit 
of M = 4 X 10^° Mq is chosen to match the z ~ 0.8 
sample that we construct below (see § 12. 2p . Approx- 
imately two-thirds of th e stars are locat ed in galaxies 
with M > 4 X 10^° Mq (|Bell et al.ll2"003[) . such that in 
terms of stellar mass our samples contain a representa- 
tive fraction of the galaxy population. 

Total magnitudes can be underestimated due to over- 
cstimatio n of the sky brig htness by the SDSS reduction 
pipeline (jLauer et al.l 1200 7). The magnitude of this ef- 
fect depends on both galaxy size and surface brightness, 
i.e., to first order on apparent magnitude. Since we are 
working with a large sample we are primarily interested 
in a statistical, magnitude-dependent correction for the 
sample as a whole. For a dozen galaxies with magni- 
tudes across the range of our sample we determine to- 
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TABLE 1 

Magnitude and Color Conversions 
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Note. — Transformations used to calculate rest-frame magnitudes and colors from 
observed magnitudes and colors for eight different redshifts. The letters v, z, and z are 
shorthand for D6O6, 2775, and 2s50, respectively. The magnitudes are total magnitudes 
for which the definitions can be found in 12.11 and 12.21 The total magnitudes go 
and Bq are not corrected for cosmologieal surface brightness dimming. Note that Bq, 
{U — -B)o, and (B — V)o are on the Vega system, whereas all other magnitudes and 
colors are on the AB system. The u — g conversion contains a correction (—0.04 mag) 
to account for the difference between the SDSS u-band zero point and the true AB 
zero point. 



tal magnitudes by han d. The true ba ckground is deter- 
mined by SExtractor (B ertin fc Arno uts 1996), using the 
global background mode, and the total flux is measured 
with the ellipse task in IRAF. We find that f/tot = 
ffmod - 0.2(15 - gmod) for gmod < 15 and gtot = .gmod for 
ffmod > 15. This corre ction is of similar m agnitude to the 
corrections applied bv lvon der Linden et a l. (2007), who 
designed a method to obtain accurate, corrected magni- 
tudes for individual galaxies. Our correction is not per- 
fect for individual galaxies but does provide a sufficiently 
accurate correction for the sample as a whole. The im- 
pact of the correction on the stellar masses that we derive 
below is ^ 0.1 dex for the brightest galaxies, such that 
even if the correction we apply is uncertain on the 50% 
level, the effect on the stellar masses is only 0.05 dex, 
which is negligible with respect to the other uncertain- 
ties on the mass estimates. The systematic error will be 
still lower. 

We adopt u-g = Umod,e - 5mod,c and g-r ^ gmod,e - 
T'mod.c as the galaxy colors, where -Umod,c, gmod,c, and 
?'mod,c are the model magnitudes, corrected for Galac- 
tic extinction. Here r,iiod,c is either the De Vaucouleurs 
model magnitude or the exponential model magnitude, 
depending on which one fits best to the observed sur- 
face brightness profile in the r band. The best-fit ra- 
dius, ellipticity, and orientation as determined in the 
r band are subsequently used to measure Umod,e and 
ffmod,e, with Only the s urface brightness as a free param- 
eter (see lLupton et al.|[2001.. for a description of the data 
reduction pipeline). The colors thus obtained are there- 
fore equivalent to the color within the r-band half-light 
radius. 

j^-corrcction s are not n egligible, even at such low red- 
shifts (Bla nton et al.l[20b3i) . We derive 5tot,o, [u — g)o, 
and {g—r)o (^-band magnitudes, and u—g and g—r colors 
if-corrected to z = 0) by followi ng the technique applie d 
to higher redshifts galaxies by iBlakeslee et al.l (|200()D . 
iHolden et al.l (jgOOG*) and Paper I. We derive linear rela- 
tions between colors and magnitudes in the observed and 
rest frames from synthet i c spec tra for stellar populations 
from lBruzual fc CharlotI ()2003D with a range of metallici- 
ties and star formation histories (exponentially declining 
star formation rates with different timescales). In Table 
[1] we show the transformations for redshifts z = 0.02, 
0.03, and 0.04. 



Subsequently, stellar masses are estimated using the 
relation between {g — r)o and M/L derived by Bel l et ahl 
((2003) for a "diet" Salpeter initial stellar mass function 
(IMF): 
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1.519(5 -r)o -0.499. 



The diet Salpeter IMF (e.g., HeireianilQQl) gives val- 
ues of M/L that are 0.15 dex lower than in the standard 
Salpeter IMF, which is due to the flat slope for stellar 
masses below 0.35 Mq. The if-corrections and color- 
M/L conversions can be combined to give a unique re- 
lation between stellar mass, apparent magnitude, color, 
and redshift: 



log 



M 

Mq 



= /(.9tot; {g - r);z) 



-0.399 



5tot -3.65(5 -r) -4.79 log 



0.03 



- 39.55 



In order to compare with previous work and the z 0.8 
sample described below, we derive the standard U — B 
and B ~ V colors and i?-band magnitudes in the Vega 
system. 

In Figure[T]we compare our stellar mass estimates with 
kinematic mass estimates, where the latter is defi ned as 
Md™/Mo = log(i?) -f 21og(crc) + 6.07 (J0rgcnscn 'et al.l 
Il996f ). Here, R is either the De Vaucouleurs or the ex- 
ponential disk scale radius in kiloparsecs as measured in 
the 5 band, depending on which profile provides the bet- 
ter fit, a nd (Te is the aperture -corrected velocity disper- 
sion fsee iBernardi et al]|2003D in kilometers per second. 
There is an average, systematic offset of 0.10 dex, the 
stellar mass estimates being lower than the kinematic 
mass estimates. This can be due to the choice of the 
diet Salpeter IMF and/or the presence of nonbaryonic 
matter, even though a is measured at the centers of the 
galaxies, which are dominated by stars. The offset varies 
slightly with mass (by < ±0.05 dex) for galaxies with 
Mphot > 4 X 10^° Mq, the mass limit of the sample used 
in this paper. Such a mass dependence of the difference 
between kinematic and stellar mass estimates is also seen 
at higher redshifts (see § 12. 2|) . The scatter in Figure [T] is 
0.17 dex, which implies that the errors in the photomet- 
ric mass estimates for individual galaxies are less than 
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Fig. 1. — Comparison between stellar masses (as used in this 
paper) and kinematic masses (oc Ra^) for galaxies at redshifts 
0.030 < z < 0.045. The dotted line indicates the one-to-one cor- 
respondence between M^y^ and A/phot- The solid line is the best 
linear fit, which has slope 0.9 and shows that the offset between 
the two mass estimates is, on average, 0.10 dex. The scatter is 0.17 
dex. 

~50%. Note that this is an upper hmit on the real un- 
certainty since kinematic mass estimates have their own 
uncertainties. We conclude that the color measurements 
and the resulting stellar mass estimates are sufficiently 
accurate for our purposes. 

The final sample consists of 2003 galaxies that are more 
massive than M = 4 x 10^° Mq. Of those, 545 are in 
the 0.035 < z < 0.045 subsample of galaxies with masses 
M > 1.25 X lO^i Mq. When we compute numbers and 
fractions of numbers of galaxies we give those high-mass 
galaxies lower weight (0.395), which is the ratio of the 
volumes occupied by the 0.035 < z < 0.045 sample and 
the 0.02 < z < 0.03 sample. 

2.2. A Mass-Selected, Volume- Limited Galaxy Sample 
in the CDF-S 

We combine photometric and spectroscopic redshift 
samples to construct a stellar mass-limited sample of 
galaxies in the CDF-S in the redshift range .6 < z < 1.0. 
We us e the photometric redshift catalog by IWuvts et al.l 
(|2007[ ). which is based on a if-band-selected sample. 
We supplemen t these with spectroscopic redshifts from 
the literature llLe Fevre et al.ll2004t iMignoH et al.ll2005l : 
Ivan der Wei et all 120051 : iVanzella et all 120061) . Only 
those galaxies for which the spectroscopic redshifts 
have been measured with reasonable confidence (best or 
second-best quality flags) are considered. Of the final 
sample, 142 out of 207 or 69% have spectroscopic red- 
shifts. We refer to the respective publications for expla- 
nations of the confidence indicators. 

Extensive comparisons by Wuyts et al. show that in 
the redshift range of our sample, there are no systematic 
differences between spectroscopic and photometric red- 
shifts, regardless of color and magnitude cuts imposed on 
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Fig. 2. — Small circles: Galaxies brighter than zgso = 24 
in the CDF-S. Large circles: Galaxies with spectroscopic redshifts 

0. 9 < z < 1.0. A mass-selected sample will be biased most strongly 
at the high-redshift end (z = 1) and against those galaxies with the 
highest -A//L, i.e., those with the reddest colors. The area indicated 
by the rectangle indicates the colors of the reddest galaxies at 2 ~ 

1. In Fig. [3] we address the redshift completeness of galaxies with 
such colors. 



the sample. There are occasional catastrophic failures, 
but their number is small (1% — 2%), and therefore this 
does not introduce significant systematic effects. 

The iC-band limit of the photometric redshift cata- 
log determines the mass completeness limit of our sam- 
ple. Since we use Vqqq-, 1775-, and Zgso-band ACS {Hub- 
ble Space Telescope Advanced Camera for Surveys) pho- 
tometry to calculate masses, luminosities, and rest-frame 
colors we quantify the completeness limit in terms of 
Z85o-band total magnitudes. The zsso-band total mag- 
nitudes are the SExtractor Best magnitude, but we add 
flux (0.2 mag) to correct for light o utside the aperture 
()Beni'tez et al.ll200l: iBlakeslee et al.l l2006V We measure 
ugoe ~ *775 a-nd 1775 — 2350 within a 0.5" diameter aper- 
ture, and we correct for differential point-spread func- 
tion (PSF) effects. The latter correction is negligible in 
^606— *775 but ranges up to 0.07 mag (for point sources) in 
*775 — -zsso- The aperture is chosen such that we measure 
the color within the effective radius or half-light radius of 
a typical z ~ 0.8 galaxy, which corresponds closely to the 
model colors we use to define the colors of the galaxies 
in the z ~ 0.03 sample (see 12. ip . 

Selection effects in a mass-limited sample will be 
strongest at the high end of a redshift bin (here 0.6 < 
z < 1.0) and for the galaxies with the highest mass- 
to-light ratio (M/L), i.e., the reddest galaxies. In Fig- 
ure [2] we show the vqqq — 1775 and 4775 — zgso colors of 
the galaxies in the CDF-S with spectroscopic redshifts 
0.9 < Zspoc < 1-0. The rectangle, which outlines the re- 
gion with the reddest galaxies at that redshift, indicates 
for which colors incompleteness first begins to affect our 
mass-limited sample. In Figure [3] we show that for those 
colors the X-band-selected photometric redshift catalog 
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Fig. 3. — Dashed line with crosses: Fraction of galaxies in the 
CDF-S with known spectroscopic and/or photometric redshifts as a 
function of zgso magnitude. The completeness begins to decrease 
at 2850 ~ 23 — 23.5. Dotted line with circles: Completeness of 
galaxies within the rectangle indicated in Fig. |2] i.e., those galaxies 
with colors similar to the reddest galaxies at z 1. This sample is 
essentially complete down to zgso = 24, which implies that we can 
use 2850 = 24 to define the mass limit of the sample. 

is virtually complete down to ^gso = 24, fainter than 
the completeness limit of galaxies with a typical color. 
In the following, we adopt zgso = 24 as the magnitude 
limit which we use to define the mass completeness limit 
of our sample. Bluer galaxies near this mass limit are 
much brighter, and incompleteness does not play a role 
for those objects. The reason for using a zgso-band mag- 
nitude limit to define the mass completeness limit instead 
of the i^-band magnitude limit of the photometric red- 
shift catalog is that we derive stellar masses by estimat- 
ing mass-to-light ratios in the rest-frame i3-band, which 
corresponds to the observed z band at 0.6 < z < 1.0. 

Stellar masses are derived with the same method as 
used for the z ^ 0.03 galaxy sample. First, observed 
magnitudes and colors (zsso, ^606 - «775, and 1775 - zsso) 
are used to derive rest-frame i?-band luminosities and 
rest-frame U — B and B — V colors. For z < 0.85 we use 
^^606 — *775 to derive [U — B)q and weoe — -z^sso to derive 
{B ~ V)q and Mb,o; for z > 0.85 we use vqog — zs5o 
to derive {U — B)q and ^775 — zsso to derive {B — V)o 
and Mb,o- These choices optimally match the observed 
and rest-frame bands, and the derived rest-frame colors 
are well behaved as a function of redshift. Examples for 
various redshifts are given in Table [TJ The rest-frame 
U-, B-, and F-band magnitudes and colors are in the 
Vega system, as opposed to the other magnitudes and 
colors, which are in the AB system. This is necessary 
to facilitate the comparison with earlier results. Finally, 
the empir i cal re lation between M/L and B — V from 
iBell et al] ()2003[ ) is used to estimate the stellar masses: 

log (^^^) = 1-737(5 - V)o - 0.994. 
For 14 early-type galaxies in the z ^ 0.8 




Fig. 4. — Rest-frame B-band luminosity vs. stellar mass for 
galaxies in the CDF-S with redshifts 0.6 < 2 < 1.0, using spectro- 
scopic redshifts if available, and, to ensure completeness, including 
galaxies with only photometric redshifts as well. Our magnitude 
completeness limit 2850 = 24 corresponds to the luminosity limit 
at 2 = 1 (dotted line). In addition, we adopt log(M/LB) = 0.8 
(dashed line) as the maximum value that the M/L of any galaxy 
can take. This is a reasonable assumption as this is the high- 
est M/L for galaxies at the high-mass, high-luminosity end of our 
sample. The magnitude and M/L limits cross at M > 4 X IQlO Mq . 
This is the mass completeness limit of our sample. At lower masses, 
our sample is incomplete for faint, red (high M/L) galaxies (see Fig. 
[Sj. The mass-selected sample contains 207 galaxies, 142 of which 
have spectroscopic redshifts. 

sample kinematic meas urements are available from 
Ivan der Wei et al] ()2005D . There is a systematic offset 
between the kinematic and stellar mass estimates of 0.1 
dex (the kinematic masses being larger) and a scatter 
of 0.3 dex. These results are consistent with earlier de- 
terminations of t he robustness of stellar mass estimates 
at z ~ 0.8 - 1 (|van der Wei et all l2006t iHolden et all 
l2006( l. Most importantly, the systematic difference be- 
tween kinematic and stellar mass estimates is the same 
(0.1 dex) for the z ~ 0.03 sample and the z ^ 0.8 sam- 
ple. Therefore, we conclude that our stellar mass esti- 
mates do not introduce systematic effects into our sam- 
ples that exceed 0.1 dex in ma ss as was found earlier b y 
both I van der Wei et al.l (|2006^ and 'HoldeneLaU (|2006f) . 
Uncertainties at that level do not significantly affect our 
results. The difference between the kinema tic and stellar 
mass estimates is mass dependent (see, e.g. jHolden et al.l 
I2OO6I and Paper I), but since this is also the case for our 
z ^ 0.03 sample (see § 12. and at the same level, this 
does not affect our results significantly. 

For late-type galaxies a comparison between kinematic 
and st ellar mass est i mates is less straightforward. Re- 
cently, iKassin et al.l ()2007| ) showed that the relation be- 
tween stellar mass and kinematic mass as defined by a 
combination of the rotation speed and the velocity dis- 
persion of the gas does not vary by more than '^0.1 dex 
over the redshift range 0.1 < z < 1.2. However, the stel- 
lar masses from Kassin et al. are determined by fitting 
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Fig. 5. — Sersic n vs. bumpiness B for all 207 galaxies in the mass-selected z ~ 0.8 sample (left) and the 2003 galaxies in the z ~ 0.03 
sample (right). The 200 galaxies with visual classifications are separately indicated, and the remainder of the sample is indicated with 
small dots. The Sersic index n is restricted to be between n = 1 and 4, but to enhance the readability of the figure we have randomly 
distributed the values of n between 0.5 and 1 for galaxies with n = 1 and between 4 < n < 5 for galaxies with n = 4. Visual classifications 
are indicated by different symbols: early types are indicated by open circles, late types by crosses. The curved line, B = 0.065(n + 0.85), 
indicates the separation between early types and late types according to the B — n classification. The correspondence between the visual 
and automated classifications is very good. See text for more details and Figure [6] for illustration. 



stellar population models to broadband spectral energy 
distributions, which is not directly comparable with our 
method. If Kassi n et al. ha d used the color-M/L re- 
lations from Bell et al.l ()2003f ). the stellar masses would 
have changed such that at 2; ~ 1 they would be ~0.2 
dex higher than those of local galaxies with the same 
kinematic mass (Paper I). We comment on the effects of 
possible biases where relevant. 

In Figure |4] we show the rest-frame i?-band luminosi- 
ties and stellar masses of the galaxies in our z ~ 0.8 
sample. The highest M/L that occurs (for luminous, red 
galaxies) is M/Lb = 6.3. We adopt this as the M/L 
upper boundary, which, combined with the magnitude 
limit derived above (2350 = 24), gives us the mass com- 
pleteness limit of our sample: M = 4 x 10^" AIq (see Fig. 
[4]). The mass- limited sample contains 207 galaxies. 

The average M/L is M/Lb = 2.4 for the galaxies 
with M > 4 X 10^° Mq. For the z - 0.03 sample 
this is M/Lb ~ 4.9. This implies 0.75 mag evolu- 
tion between the two samples in the B-band, or 1 mag 
per unit redshif t. This is consist ent with the evolu- 
tion of L* (e.g.. iBrown et al]l2007f l. but somewhat less 
than the M/Lb evolution inferred via fundamental plane 
analyses for z 1 field early-type galaxies with simi- 
lar rnasses (0.55 dex: Ivan der Wei et al.ll2005l: lTreu et all 
120051 : Ivan Dokkum fc van der Mare]|l2007f ). This differ- 
ence may contain contributions from a systematic uncer- 
tainty in the high-z mass estimates as described above, 
and an increased number of very red, dusty galaxies at 
z - 1 with high M/L. 

2.3. Morphological Classifications 

Traditionally, visual morphological classifications are 
used for morphological studies of high-rcdshift galaxy 
populations, for example the evolution of the MDR 



(e.g.. iPostman et"al] l2005f ). In recent years, however, 
increased sample sizes have initiated the emergence 
of several alternative, automated classification schemes 
(■Conselice et al.. 2000 : Abjraham et al. 2003; Lotz et al] 
12004 iBlakeslee et al.H2006[ ). In addition to the tradi- 
tional, visual classifications we therefore also deploy a 
quantitative measure of the galaxy morphologies. This 
enables us to analyze the large sample of z ~ 0.03 galax- 
ies and to examine, in a quantitative manner, the effect 
of cosmological surface brightness dimming on the de- 
tectability of disks, spiral arms, and irregularities. The 
visual classifications are used to verify the robustness 
of this automated classification method. In this pa- 
per we simply distinguish between early- and late-type 
galaxies, or, in the T-type classification scheme (see 
IPostman et al.ll2005f l. between galaxies with type T < 
(E, E/SO, and SO) and galaxies with type T > (Sa and 
later). The separation into subclasses, for example E and 
SO, will be investigated in a future paper. 

Visual morphologies were determined by A. v. d. W. 
by examining the g-band images of 200 randomly chosen 
galaxies in the mass-selected z ^ 0.03 sample and the 
zgso-band images of all 207 galaxies in the z ~ 0.8 sam- 
ple. For 35 galaxies in the z ~ 0.8 sample with T 
or peculiar morphologies M. P. also examined the im- 
ages and assigned morphologies. This provides an es- 
timate of the random uncertainty and e nsures that the 
same threshold (T = 0) is adopted as in IPostman et al] 
()2005f ). Based on these results, we conservatively as- 
sume a 5.5% uncertainty in the early-type galaxy fraction 
due to classification errors. We note that the Z85o-band 
imaging with an integration time of 12,000 s and a IOct 
point-source detection limit of zsso = 27 is very deep 
with respect to the faintest galaxies in our sample with 
Z850 = 24. In Figure [5] we indicate early-type galaxies 
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with open circles and late-type galaxies with crosses. 

Quantitative morphologies are dete rmined for all 
galaxies with the technique outlined by iBlakeslee et al.l 
|^0(f), who showed that the Sersic parameter n 
and bumpiness parameter B can effectively distin- 
guish between ea rly- and late-type galaxies. GALFIT 
(|Peng et"al1l2002D is used to fit Sersic profiles to the g- 
band and zsso-band images; n is constrained to values 
between n = f and 4. The only reason for not allowing 
larger values of n is that the effective radii for high values 
of n become large and uncertain. This would strongly af- 
fect the measurement of the bumpiness B, which is the 
rms of the residual as measured within two effective radii. 
In order to reduce the effect of shot noise the residual is 
slightly smoothed. In addition, other objects and a circu- 
lar region around the center are masked. This central re- 
gion is masked out because of slight PSF mismatches and 
because of central, small-scale deviations from a Sersic 
profile. Because of the latter effect the adopted radius is 
distance dependent. On visual inspection we choose a 4 
pixel radius for galaxies at z = 0.03 in the low-z sample 
and a 2 pixel radius for galaxies aX z = 0.80 in the high-z 
sample, both scaled inversely proportional with angular 
diameter distance for galaxies at different redshifts. 

The advantage of the B — n classification over other 
automated classifiers is that PSF smearing is taken into 
account. This is essential, in particular when redshift- 
dependent trends are investigated in data sets with very 
different photometric properties. In Figure[6]we show im- 
age cutouts and give the n and B parameters for 20 ran- 
domly selected galaxies from our z ~ 0.0 and 0.8 samples, 
in order to illustrate the correspondence between the au- 
tomated classification method and their morphological 
appearance. The parameters n and B are shown for all 
galaxies with visual morphologies in Figure [51 A sim- 
ple linear relation separates early- and late-type galaxies: 
B = 0.065(n + 0.85) (Fig. El solid line). No t e that these 
coefficients are taken froni IBlakeslee et al.l ()2006t ). and 
are not chosen to optimize the correspondence between 
the visual and B — n classifications for these particu- 
lar samples. The correspondence between our results, 
in particular for the z ~ 0.03 sample, and the results 
by Blakeslee et al. illustrates the power of the B — n 
method. Once PSF smearing effects and noise properties 
are properly taken into account, data sets with roughly 
similar spatial resolution provide a unique measurement 
of galaxy morphology through the B — n classification. 
This extends the conclusions from Blakeslee et al. who 
could not claim universality of their selection criteria be- 
cause of their homogeneous data sets and the small red- 
shift range (z = 0.83 — 0.84) of their samples. 

The visual and B — n classifications agree for 93% of 
the galaxies in the z ^ 0.8 sample (193 of 207) and for 
89% of the galaxies in the z ~ 0.03 sample (179 of 200; 
see also Fig. |S]). Half of the mis-classifications have n- 
and i3-values that put them close to the criterion that 
separates the early types from the late types (Fig. |S1 
solid line). The disagreement can be explained by the 
Hmited signal-to- noise ratio [S/N) of the images, which 
shows that the true random uncertainty in the visual and 
B ~ n methods is ~4%. The other half of the disagree- 
ments are mostly due to strong, central deviations from 
a Sersic profile (probably point sources) or, for several 
galaxies in the z ^ 0.03 sample, large-scale deviations 



from a smooth profile that our fitting method does not 
place within two effective radii. Despite the nonnegligi- 
ble numbers of erroneous classifications (~10%), the net 
difference between the two classification methods in the 
ratio of the numbers of early-type and late-type galaxies 
is less than 0.5%. This is true for both the z ~ 0.03 
sample and the z ~ 0.8 sample. Hence, there is virtually 
no systematic difference between the visual and B — n 
classifications for the samples as a whole. 

Even though the z ~ 0.03 and 0.8 samples have in- 
ternally consistent morphological classifications, there 
may be a systematic difference between the two sam- 
ples, which are located at very different cosmological dis- 
tances. Our morphological classifications could depend 
on redshift, mainly because of the lower S/N of the high- 
z galaxy images. Because the scales of the PSF and 
pixels of the SDSS and ACS data sets are very similar in 
terms of physical size at the sample redshifts, we can test 
the redshift dependence of our morphological classifica- 
tion. We take 50 galaxies in our sample spanning the full 
range in magnitude and with redshifts 0.025 < z < 0.030, 
such that the SDSS pixel scale corresponds to the pixel 
scale of an ACS image of a galaxy at z = 0.80. We add 
noise to the 5-band SDSS images, taking into account 
i^-corrections and the relative depths of the SDSS and 
ACS images, and reapply our two classification meth- 
ods, which gives no. so and -Bq.sOj the values for n and B 
as inferred from the simulated images. In Figure [71 we 
compare those with the original values of n and B, and 
we find that neither quantity shows a systematic differ- 
ence. In addition, the visual morphologies of the sim- 
ulated images are not different from the original visual 
morphologies: all features that decide on the morphology 
of the z 0.03 galaxies are still visible in the simulated 
z = 0.80 images. We conclude that our classification 
methods do not suffer from systematic effects that intro- 
duce significant differences between the morphologies of 
the z ~ 0.03 and z 0.8 galaxies. 

2.4. Local Density Estimates 

Traditionally, local projected surface densities are com- 
puted by measuring the distance to the n-th nearest 
neig hbor brighter than, for example, Mv-l-0.8z < —19.78 
(e.g. JPostman et aLll2005t ). However, since we work with 
mass-selected samples in this paper, it is more consistent 
to measure the distance to the 71th nearest neighbor that 
is more massive than, in this case, M > 4 x 10^° Mq. Be- 
low we discuss the difference between the two approaches. 

We derive the local galaxy densities for the 142 galax- 
ies in our z ~ 0.8 sample with spectroscopic redshifts. 
We compute the distance to the n-th nearest neighbor 
with a difference in rest-frame radial velocity of Av < 
1000 km s~^ and more massive than M = 4 x 10^° Mq 
(or, alternatively, more luminous than My + 0.8z < 
— 19.78). Here n is the number of galaxies with such 
properties in the CDF-S, with a maximum of 7. The 
reason for allowing n to be smaller than 7 is the limited 
number of galaxies at a given, spectroscopic redshift. We 
note that for galaxies with fewer than 7 neighbors (12 
of 142 for the luminosity-limited sample, 36 of 142 for 
the mass-limited sample) the density estimate does not 
systematically change for any value 3 < n < 7. The in- 
ferred surface densities are corrected for edge effects and 
for incompleteness of the spectroscopic redshift catalog. 
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Fig. 6. — Top two rows: 3.4" X 3.4" ^gso-band cutouts for 10 randomly chosen galaxies from the sample at 0.6 < z < 1.0. Bottom 
two rows: 51" X 51" g-band cutouts for 10 randomly chosen galaxies from the sample at z ~ 0.03. The physical scale of both sets of 
images is comparable: ~30 kpc. Visual classifications (late or early) and n and B parameters are listed for each galaxy to illustrate the 
correspondence between the visual morphologies and the automated classification method (see also Fig. [Sj- 



as inferred from the spectroscopic + photometric rcdshift 
catalog. 

The distributions of Sm and El are shown in Figure 
[HI Sm is typicahy half of El , which is only a small differ- 
ence considering the large spread (3 orders of magnitude) 
of densities in the samples. For galaxies with only few 
neighbors, the computed values are likely upper limits 
due to small area of the CDF-S. Obviously, for galaxies 
without close neighbors, the computed densities are cer- 
tainly upper limits. The densities as computed here are, 
in principle, not directly comparable wit h the computa- 
tions for galaxies in cluster environments (jPostman et al.l 
[200 5) as those use a wider rcdshift range to search for 
neighbors to compensate for the large velocity disper- 
sion of the cluster members. A systematic difference of 
up to a factor of 2 may be expected due to the difference 
in bin width. However, since the environments under 
consideration have a range in density of several orders of 



magnitude this does not compromise our inferences. 

In a method analogous to that used for the z 0.8 
sample, we also compute local galaxy densities by find- 
ing the distance to the seventh nearest neighbor with 
M = 4 X 10^° Mq (or My + 0.8z < -19.78) and 
Av < 1000 km for each galaxy in the z 0.03 sam- 
ple. The spectroscopic survey is 100% complete at these 
magnitudes, therefore corrections are not needed. The 
distribution of E is shown in Figure [H As is the case 
for the z 0.8 sample, Em and El typically differ by a 
factor of 2 for the z ~ 0.03 sample (see Table [1]). 

There is a difference of ^0.7 dex between the me- 
dian values of E for the z ^ 0.03 and 0.8 samples. 
This is mainly because we use the angular diameter dis- 
tance Da to calculate the physical distances between 
the galaxies, and not the transverse co-moving distance 
Dm = Da{1 + z). This provides, to first order, rcdshift- 
independent values of the local density for bound systems 
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Fig. 7. — Sersic parameter n and bumpiness parameter B for 
50 galaxies at rcdshifts 0.025 < z < 0.030 as determined from the 
original SDSS g-band images vs. the difference between n and B 
as determined from simulated z = 0.80 z-band images of the same 
galaxies and the original values. There are no systematic differ- 
ences between the original low-z and simulated high-z morphologi- 
cal indicators, which indicates that the morphologies as determined 
for our 2 ~ 0.03 and ~ 0.8 samples are not systematically different. 

such as clusters, but for expansion-dominated regions the 
density will be redshift dependent. Using Dm would pro- 
vide a redshift-independent measurement of the density 
and would lower the z 0.8 densities by 0.5 dex with 
respect to the z ~ 0.03 densities. This demonstrates 
that the galaxies in both samples are located in similar 
(co-moving) environments, besides an average, apparent 
over-density of a factor of 1. 5 — 2 in the CDF-S which 
has been noted before (e.g.. iBundv et alll2005[ ). This 
difference is not of great interest for this study since we 
consider a range in density that spans several orders of 
magnitude (see § 13. 2p . Since we are interested in the 
properties of the galaxy population as a function of phys- 
ical density, and to facilitate the comparison with work 
on cluster galaxies, we use E as calculated with Da in 
the remainder of this paper. 

3. THE EVOLUTION OF THE EARLY-TYPE GALAXY 
FRACTION AND THE MORPHOLOGY-DENSITY 
RELATION 

In § 13.11 we calculate the evolution of the early-type 
galaxy fraction in our field samples. In § 13.21 we com- 
pare this with samples of cluster galaxies over the same 
redshift range and determine the evolution of the MDR. 

3.1. The Field Early-Type Galaxy Fractions at z 0.8 

and 

Our z ~ 0.03 sample consists of 2003 galaxies, the 
z ~ 0.8 sample of 207 galaxies. These volume-limited 
samples are complete down to a stellar mass limit of M = 
4 X 10^" Mq, about half of the typical mass of a local field 
galaxy, or Mb ^ —20.1 for red galaxies in the z ^ 0.8 
sample. This luminosity limit for red galaxies is similar 
to the z = 0.8 - 1 luminosity limits of the COMBO- 




FlG. 8. — Distribution of local surface densities of the galaxies in 
our samples, (a) Distribution of the 0.6 < z < 1.0 sample, where 
the vertical dashed line indicates the value of E corresponding to 
one galaxy within the area covered by the survey data: objects 
with values of E close to this line should be regarded as upper 
limits, (b) Same as (a), but for the 0.015 < z < 0.045 sample. The 
solid histograms show the distributions of Sm , the surface density 
as determined by regarding neighbors that are more massive than 
M = 4 X 10^" Mq . The dotted histograms show the distributions 
of El, the surface density as determined by regarding neighbors 
that are brighter than My + 0.8z = -19.78. 



17 (iBeh et al.ll2004hD. pEEP2 (|Faber et al.l [200l and 
NDWFS (jBrown et al.l I2007D data sets used for recent 
red-sequence studies. 

The early-type galaxy fraction at z 0.8, for galax- 
ies that are more massive than M — 4 x 10^° Mq, is 
/et — 0.48 ± 0.07. This is based on the quantitative 
morphological classification method described above, but 
we find the same if we use the visual classifications (see 
also ij 12. 3p . If only galaxies with spectroscopic redshifts 
are included (142 of 207 have spectroscopic redshifts), 
we find /et = 0.52. The reason that the spectroscopic 
sample has a slightly, but not significantly, higher early- 
type galaxy fraction is that several of the spectroscopic 
campaigns specifically targeted red objects. 

The main contributor to the error of 0.07 is the uncer- 
tainty in the morphological classifications (0.055). The 
second main source of uncertainty is Poisson noise due 
to the limited sample size (0.036). Smaller contributions 
include errors in photometric redshifts, color transforma- 
tions, and M/L estimation. All these individual contri- 
butions are added in quadrature to obtain 0.07. 

The early-type fraction at z ~ 0.03, for galaxies that 
are more massive than M = 4 x 10""^° Mq, is /et = 
0.43 ± 0.03. Essentially, the only significant contributor 
to the error is the uncertainty in the morphological clas- 
sifications, which is slightly smaller than for the z ^ 0.8 
sample because of the higher S/N of the images. 

The early-type galaxy fractions in the mass-selected 
z ~ 0.03 and 0.8 samples are not significantly different 
(see Table [T] and Figure [5]). This implies that if only 
two morphological classes are considered, then the mor- 
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TABLE 2 
Sample Characteristics 







/et 


/rs 


/n>2.5 


Med{SM) 


Med(SL) 




< SFR> 


Med{SFR) 


z 


N 


(%) 


(%) 


(%) 






Med(l/ - B) 


M© yr 


Mq yr 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


0.020 <z < 0.045 


2003 


43 ±3 


68 ± 1 


64 ±2 


0.17 


0.4 


0.40 


3.1 ± 1.5 


0.4 ±0.2 


0.6 < 2 < 1.0 


207 


48 ± 7 


64 ±3 


72 ±4 


1.0 


2.1 


0.23 


12 ±6 


<3 



Note. — Col. (1) Redshift. Col. (2) Number of galaxies in the sample. Col (3) Early-type galaxy fraction as inferred 
from the automated B — n classifications (see ^ |2.3l l . Col. (4) Fraction of galaxies located on the red sequence as defined in 
^ 14.21 Col. (5) Fraction of galaxies with Sersic n parameter larger than n — 2.5 I2.3| l . Col. (6) Median local surface density 
of galaxies with masses M > 4 X 10^" Mq (scg §[23J. Col. (7) Mod ian local surface density of galaxies with luminosities 
Mv + 0.82 < -19.78 (see §[2T3}. Col. (8) Median rest-frame U - B color as computed in ^ g [27T1 and [2T2I Col. (9) Moan 
SFR (§ I4.lt : (9) Median SFR (§ I4.lt . The errors on the SFRs are assumed to be 50%, mainly to account for systematic 
uncertainties. 



phological composition of the galaxy population did not 
change significantly between z ~ 0.8 and the present, at 
least above our mass limit of M = 4 x 10^° M©. If any- 
thing, the early-type fraction is higher at z ~ 0.8 than 
at z ~ 0.03. 

3.2. Comparison with Cluster Calaxies: The 
Morphology-Density Relation 

As we described above, the early-type galaxy fraction 
in our field samples is 40% — 50% at redshifts < z < 1. 
In Figure ini we compare this with the results from Paper 
I, where we measure the evolution of the cluster early- 
type fraction for galaxies down to the same mass limit 
(M = 4 X 10^" Mq) and over the same redshift range. 
For galaxies more massive than M = 4 x 10^° Af© the 
early- type galaxy fraction is constant with redshift, both 
in the field and in massive clusters (see Fig. [9]). The same 
MDR exists at z = and 0.8, with early-type fractions 
of 40% — 50% in low-density environments and >80% in 
high-density environments. 

The cluster samples from Paper I consist of galaxies 
with local densities > 50 Mpc"^, our field samples consist 
of galaxies with densities < 50 Mpc~^. In Figure [TOl we 
show the MDR at redshifts z ~ 0.03 and - 0.8. In 
this figure we split our z ~ 0.03 field sample of 2003 
galaxies into five density bins, ranging from 0.01 to 10 
Mpc~^. Our z ~ 0.8 field sample is represented by a 
single data point, because of the limited sample size. The 
average logarithm of the density of the 142 galaxies with 
spectroscopic redshifts is used to estimate the density 
for the sample as a whole. The MDR has not evolved 
significantly since z 0.8 for galaxies more massive than 
M = 4 X lO^*' Mq. The early-type galaxy fraction in 
regions with a given local density has remained constant 
over the past ~7 Gyr. 

At first sight, this result seems to be at odds 
with previous studies that reported significant evolu- 
tion of the early-typ e galaxy fraction (jSmith et al.ll2005t 
iPostman et al.ll2005l ). However, those analyses are based 
on luminosity-selected samples, and, as is shown in Paper 
I, the observed evolution is driven by blue, star-forming, 
low-mass galaxies. Samples that are complete in terms 
of stellar mass do not show evidence for evolution in 
the early-type fraction. The evolution in the early-type 
fraction seen in luminosity-selected samples is driven by 
the evolution of the fraction of SO galaxie s, whereas the 
fraction of E galaxies is roughly constant ("Dressie r et al " 



1997t iTreu et al.ir2003l : IPostman et"an[2005. : ,Desai et al 



E and SO galaxies evolve in mass-selected samples. 

The MDR in the z 0.03 data set shown in Figure [TOl 
suggests that the early-type fraction is low (~40%) for 
all local surface densities Em ^ 1 Mpc~^. The transition 
to very high early-type fractions (~100%) happens at 
1 Mpc"^ < Sm < 100 Mpc"^. This suggests that there 
is a critical density above which galaxies are progres- 
sively more likely to undergo morphol ogical transforma- 
tions, as was earlier noted bv lGoto"er al. (2003). Despite 
the small area and the small dynamic range in densities, 
there is evidence for this transition in our z ~ 0.8 sam- 
ple: if we take the galaxies with spectroscopic redshifts 
in the highest 20th percentile of the density distribu- 
tion (typically, Sm ~ 30 Mpc~^), we find an early- type 
fraction of /et = 0.8 ± 0.1; for galaxies in the lowest 
20th percentile (typically, Sm ~ 0.2 Mpc"^), we find 
/et = 0.4 ±0.1. Without over-emphasizing these num- 
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Fig. 9. — Early-type galaxy fraction as a function of redshift. 
The circles indicate our samples of field galaxies (S < 50 Mpc~^) 
at z ~ 0.03 {open blue circle) and 2 ~ 0.8 (filled red circle). The 
squares indicate the samples of cluster galaxies (S > 50 Mpc~^) 
from Paper I at z = 0.02 (open blue square, representing the Coma 
Cluster) and z = 0.83 (filled red square, representing the CL 0152 
and MS 1054 clusters). The fraction of early- type galaxies has not 
changed significantly since z ~ 0.8 for galaxy populations in similar 
environments. 
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Fig. 10. — Morphology-density relation for mass-selected galaxies (M > 4 X 10^" Mq) at z ^ and z ^ 0.8. The symbols are the 
same as in Figure [9] Our z ~ 0.03 field sample (open blue circles connected by the dotted line) have been split into five density bins. 
Over almost 3 orders of magnitude in density the morphology-density relation has been in place since at least z = 0.8 and has not evolved 
within the measurement errors. The difference between the two panels is the use of Sm (the surface density of galaxies more massive than 
M = 4 X IQl" Mq, or ~ 0.5M.) in (a), and the use of Sl (the surface density of galaxies more luminous than My + O.Sz < —19.78, or 
~ 0.5Lt) in (ft). Other than a shift of 0.3 dex in density the use of Sm and El is interchangeable. 



bers because the sub-sample with spectroscopic redshifts 
is not necessarily representative of the whole z ~ 0.8 
sample, we do note that these values follow the same 
trend as seen in the z ~ 0.03 sample (see Fig. [TO)) . 

4. DISCUSSION 

4.1. The Evolution of the Star Formation Rate 

Our conclusion that for galaxies with masses M > 
4 X 10^*^ Mq the early-type galaxy fraction has not signif- 
icantly changed in any environment over the past 7 Gyr 
has to be reconciled with the well-constrained decrease 
of the star formation rate (SFR) d ensity by an order 
of magnitude since z ^ 1 (s ec, e.g., 'Madau et al."1996l: 
Wolf et al. 2003; Le Floc'h et al. 2005; Bell ct al. 2005; 
Perez-Gonzalez etirn2005t IZheng et al.l i2006'). It may 



seem counter-intuitive to have the SFR decrease rapidly 
without a change in the morphological mix of the galaxy 
population. On the other hand, a non-changing early- 
type galaxy fraction does not necessarily imply that the 
galaxy population itself does not evolve. Below we inves- 
tigate the star formation activity in our samples. 

We use the very deep, publicly available GOODS^ 
24 /im imaging from t he Multiband Imaging Photome- 
ter for Spitzer (MIPS; iRieke et all [2004) to verify the 
star formation activity in our z 0.8 sample of mass- 
selected galaxies. Total 24 fj,m fluxes are measured by 
fitting the MIPS PSF to the positions of objects detected 
in the 3.6 Atm ban d of the Infrared Array Camera (IRAC; 
iFazio et al.l 12004} ) on board Spitzer. This technique re- 
duces the effect of blending and allows for more accu- 
rate photometric measurements (for more details, see 
iWuvts et al.ll200l . The 24 /im flux is used to estimate 

* The Great Observatories Origins Sur- 
vey' (GOODS) data are available at 

|http://data.spitzer. caltech.edu/popular/goods/Documents/goods-dr3, 



the bolom etric infrared luminosi ty Xir via the conver- 
sions from lCharv fc ElbazI (l200ll). Subseque ntly, we use 
LiR to estimate the SFR ( KennicuttI 119981 ). where the 
SFR is scaled down by 0.15 dex to account for the dif- 
ference between the standard Salpeter IMF and the diet 
Salpeter IMF that we use. For more details concerning 
the iiT-correction and unc ertainties (roughly a factor of 
2 in the SFR) we refer to Ivan der Wei et all (|2007. ). We 
conservatively use 50 /uJy as the flux limit, as at lower 
flux levels the poor constraint on the positions of 24 fim 
sources requires verification of the optical counterpart for 
each individual object, which is beyond the scope of this 
paper. 

In Figure [TT] we show the SFR and the SFR per unit 
stellar mass, the specific SFR, as a function of rest-frame 
U - B color for the 187 (of 207) z ~ 0.8 galaxies with 
MIPS coverage. Of the early-type galaxies, 11% are de- 
tected at 50 /jJ y, at least half of which are likely active 
galactic nuclei (Ivan der Wei et aLllMfTl ). Of the late- 
type galaxies, 76% have fluxes > 50 /.tJy, the majority of 
which are most likely due to star formation. The typical 
SFR of a late-type galaxy is ~ 20 Mq yr~^ , whereas that 
of an early type is < 5 AIq yr~^. The specific SFRs are 
< 10~^° yr~^ for early types and < 10~^ yr~^ for late 
types. These values imply an increase in stellar mass of 
less than 10% per Gyr for the early-type galaxies and 
up to 100% for the late types (see Fig. [Tl]). These find- 
ings are consistent with the work by Bell et al. (2005), 
who show that the star formation activity at z 0.7 is 
mainly due to star formation in normal spiral galaxies 
w ith masses > 2 x 1 ^° Mq . 

'Brinchm ann et al] (|2004f ) derived SFRs for galaxies in 
the SDSS by modeling their optical emission lines. The 
most robust determination of the SFR for nearby galaxies 
such as those in our z ~ 0.03 sample is that measured 

html| 
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Fig. 11. — Top left: SFR (derived from 24 /im photometry) as a function of rest-frame color for the z ~ 0.8 sample. The red circles are 
early-type galaxies, the blue circles are late-type galaxies. Filled circles are objects with significant (> 5a) 24 fim detections, open circles 
are objects without significant detections (upper limits are shown). Bottom left: For the same sample as in the top left panel, the specific 
SFR as a function of color. Top right: SFR (derived from modeling optical emission lines; Brinchmann et al. 2004) as a function of color 
for the z ~ 0.03 sample. Open red circles are early-type galaxies, and filled blue circles are late-type galaxies. Bottom right: For the same 
sample as in the top right panel, the specific SFR as a function of color. For both the low- and high-redshift samples, the star formation 
activity is low in early-type galaxies. SFRs are higher for late types, with a marked decrease in SFR for red late-type galaxies between 
z = 0.8 and the present. The dotted lines in the two bottom panels indicate the expected specific SFRs associated with galaxies that will 
double/have doubled their stellar masses between 2 = 1 and the present under the assumption that their SFRs decline (exponentially as 
a function of redshift) by an order of magnitude between 2 = 1 and the present. The decreased star formation activity in our z ^ 0.03 
sample with respect to our 2 ~ 0.8 sample is consistent with the decline of the cosmic SFR. 



within the spectroscopic fiber of the SDSS spectrograph. 
We use t he stellar mass derived within the same fiber 
aperture (iKauffmann et all l2003f ) and scale the SFR, 
comparing the fiber stellar mass with our total stellar 
mass (§ 12. ip . We show the SFR and specific SFR in 
Figure [TTl The global trend is similar to that observed 
at z 0.8, with early-type galaxies showing lower star 
formation activity than late- type galaxies. 

It should be kept in mind that SFRs of the z ^ 0.03 
galaxies are determined from optical emission lines as 
measured with a 3" fiber, such that extended star forma- 
tion can be missed, and a systematic difference between 
the total SFRs of the z ~ 0.8 galaxies is created. In addi- 
tion, it is well known that IR-derived and emission- line- 
derived SFRs can be intrinsically different and depend 
on metallicity and the amount of extinction. Therefore, 
we have to restrict ourselves to the description and in- 
terpretation of evolutionary trends of at least an order 
of magnitude. 

Despite the large uncertainties and the high upper limit 
of 5 Mq yr~^ for the z ^ 0.8 sample. Figure [Til shows 
that the overall star formation activity decreases between 
z ^ 0.8 and the present (see also Table [Ij. One promi- 
nent difference between the z ~ 0.03 and ^ 0.8 samples 
is the marked decrease of star formation activity by 2 or- 
ders of magnitude in late-type galaxies that are located 
on the red sequence. In the local universe, their SFRs are 
as low as those of early-type galaxies, whereas a.t z 0.8 
their SFRs are at least an order of magnitude higher than 

^ The stellar masses and SFRs are publ icly available at 
|http://www. mpa-garching.mpg.de/SDSS/DR4/ 1 



those of the early types. 

In Figure [Tl] we indicate the specific SFR expected 
for galaxies that double their stellar mass between z = 1 
and the present, assuming an exponential decline in SFR 
by a factor of 10 over that period. Such a level of star 
formation, both at z ~ 0.03 and at z ~ 0.8, coincides 
with the typical specific SFR of late- type galaxies. The 
decrease in star formation activity as observed in our 
samples follows the cosmic average trend of decline by 
an order of magnitude since z ~ 1. At the same time, 
the morphological mix has remained unchanged, at least, 
if only two morphological classes are considered. These 
results are therefore not mutually exclusive and have to 
be reconciled with each other. 

4.2. The Red Sequence and the Evolution of the Stellar 
Mass Density 

In Figs. [12] we show the rest- frame U ~ B color distri- 
bution as a function of stellar mass for the low- and high- 
redshift samples. The U — B colors of the z ~ 0.8 sam- 
ple are corrected for evolution within the sample redshift 
range 0.6 < z < 1.0, assuming A(f7 ~ B) ^ 0.22(z - 0.8) 
(see below). This is a small effect of typically 0.02 mag. 
We define the early-type galaxy red sequence for the 
z ~ 0.03 sample by a linear fit, with the slope and the 
zero point as free parameters, and iteratively rejecting 
3cr outliers (8% are rejected, virtually all blueward of the 
red sequence). In order to minimize systematic effects, 
we define the red sequence for the z ~ 0.8 sample with 
the same slope (0.08 mag dex~^) and scatter (0.05 mag), 
with only the zero point as a free parameter. The red 
sequence for both samples is shown in Figure [T2l Color 



VAN DER WEL ET AL. 



13 



0.6 



0.4 



£0.2 




03 

I 



-0.2 - 



-0.4 



T 1 — I — I — r-r 



0.6 



0.4 



D On " • o o* ^„ • 



0.2 



-0.2 - 



-0.4 



-1 — I — I — r-r 



10' 



10' 



10' 



10' 



Fig. 12. — Stellar mass vs. rest-frame U — B color for the z ~ 0.8 sample (ie/t) the z ~ 0.03 sample {right). Early-type galaxies are 
indicated by filled red circles symbols, late-type galaxies by open blue circles. Color evolution is readily visible, with the red sequence being 
significantly bluer at 2 ~ 0.8 than at 2 ~ 0.03. The solid lines indicate the best fits to the red sequence as defined by the early-type galaxies, 
and the dashed lines separate blue and red galaxies (see text for an explanation). The relative numbers of red galaxies in the z 0.03 
and ^ 0.8 samples are not significantly different (see also Table [TJ- In addition, the fraction of red galaxies with late-type morphologies is 
similar in both samples. 

tirely appropriate to compare galaxy samples at different 
redshifts down to the same mass limit. If galaxies evolve 
by a factor of 2 in stellar mass between z = 0.08 and 
the present, the mass limit ai z — 0.8 should be cho- 
sen 0.3 dex lower than at z = 0, otherwise many of the 
progenitors in the z — Q sample will not be included 
in the = 0.8 sample. We test this effect by choosing a 
mass limit for our z ~ 0.03 sample that is 0.3 dex higher 
(M = 8 X lO^" Mq) than for the z - 0.8 sample. The 
early-type fraction increases slightly, from 43% ± 3% to 
48% ± 3%. Due to the errors it is not clear whether this 
agrees better with the z ^ 0.8 early-type galaxy fraction 
(48% ± 7%), but it is striking that the values are the 
same. 

We conclude that the growth of the red-sequence 
galaxy population must happen in such a way that the 
MDR is conserved. The accretion of red galaxies onto 
the red sequence and/or the growth of galaxies on the 
red sequence must preserve the relative number of late- 
and early- type galaxies with masses M > 4 x 10^" Mq. 
The galaxy population evolves along the MDR, whereas 
the MDR itself is a fixed relationship between environ- 
ment and morphological appearance. 

Since many star-forming, late- type galaxies are located 
on the red sequence, the stellar mass density of red- 
sequence galaxies will increase through in situ star for- 
mation. The question is how much this contributes to the 
observed increase in stellar mass density of red galaxies. 
If we make the assumption that the SFR of the galax- 
ies in our z ~ 0.8 sample decreases exponentially at the 
same rate as the cosmic average, i.e., by an order of mag- 
nitude between z = 1 and the present, the average late- 
type galaxy on the red sequence will increase its stellar 
mass by 20% — 80%. The large uncertainty is caused 
by the absolute uncertainty in the IR-derived SFR of a 
factor of 2 (see § 14. ip . Given such an increase in stellar 



evolution between z ~ 0.8 and ~ 0.03 of 0.15 mag, or 
0.20 mag per unit redshift, is readily visible. We also fit 
the slope and scatter of the z ^ 0.8 red sequence inde- 
pendently, resulting in a higher scatter (0.08 mag) and a 
steeper slope (0.11 mag dex~^). 

For both samples we define blue galaxies as galaxies 
with colors more than 2 a below the red sequence (with 
fj = 0.05 as defined above; see Fig. [T^ . With this defini- 
tion, wc find similar red-galaxy fractions of 64% and 68% 
for the z ^ 0.8 and ~ 0.03 samples, respectively (see Ta- 
ble [T|). Like the early- type fraction, the red-galaxy frac- 
tion shows no strong evolution with redshift. The red- 
galaxy fraction is significantly highe r than the early-typ e 
fraction in both samples (see also iBundv et al.l (20061 ) . 
This is due to the presence of many late-type galaxies on 
the red sequence: 38% ± 8% for the z ^ 0.8 sample and 
45% ± 3% for the z ~ 0.03 sample of the red galaxies 
are late types (see also Fig. [13]). Conversely, in either 
sample less than 20% of the blue galaxies have early-type 
morphologies. 

Given the fact that the majority of the galaxies in 
both our samples are located on the red sequence, our 
earlier conclusion that the early-type galaxy fraction at 
a given density does not evolve with redshift seems, at 
first sight, to be at odds with the observation that the 
stellar mass densit y of red-sequence galaxies has d oubled 
since z = 1 (e.g.. [Bell et"al]|2004bt iBrown et al.ll2007t ). 
The luminosity limit for red galaxies at z = 1 in our 
sample (Mb ^ —20) is s imilar as that of, e.g.. lBell et al.l 
(|2004b[) and IBrown eTal . (2007); therefore, the increase 
in mass density found by those authors has to take place 
above the mass threshold of our samples. In fact, even in 
our samples, which cover only small volumes, we observe 
an increase in stellar mass density from z = 0.8 to the 
present. 

If the mass function of galaxies evolves, it is not en- 
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Fig. 13. — Rest-frame U — B color distribution of the mass selected z ~ 0.8 [left) and z ~ 0.03 (right) samples. The solid lines indicate 
the early-type galaxies, and the dashed lines indicate the late-type galaxies. The relative numbers of red galaxies and the relative numbers 
of red galaxies with late-type morphologies are not significantly different at 2 ~ 0.03 and ~ 0.8. 

mass, there will be galaxies that are below our mass cut- 
off at z 0.8 but would not be during the present epoch. 
We can estimate this effect by calculating the total stel- 
lar mass in red, late-type galaxies in our z ~ 0.03 sample 
with masses less than 1.2 — 1.8 times the threshold. Their 
contribution to the total stellar mass in red galaxies is 
2% — 11%, such that the total increase in the stellar mass 
density through in situ star formation in red late types 
is a factor of 1.2 — 2. This is close the observed evolution 
of a factor of ^2. One caveat is the possible evolution 
with redshift of the color-M/L relation, which we use to 
estimate the stellar masses (see § 12. 2[) . This may cancel 
out the growth in stellar mass as estimated above. 

With the hypothesis that merging does not play a role 
in shaping the red sequence galaxy population, the above 
described process of star formation in late types, but 
not in early types, will change the morphological mix of 
a mass-selected sample of red galaxies. The early-type 
fraction will then decrease with time or increase with 
redshift. As mentioned in § 14.21 the fraction of early- 
type galaxies on the red sequence is indeed consistent 
with an increase from z ~ 0.03 to ^ 0.8 (from 55% to 
62%). Even though this difference is not significant, at 
least it is consistent with a growth in stellar mass of late- 
type galaxies by a factor of 1.5. 

These speculations depend on the assumption that a 
large reservoir of lower mass, red galaxies is available 
to provide the necessary growth of stellar mass through 
star formation on the red sequence. Given the evi- 
dence for the opposite, i.e., the lack of faint, red galaxies 
at significant look-ba ck times (e.g.. I^naka et ani2005l : 
iDe Lucia et alj [20071 ). truncation of star formation in 
blue galaxies is possibly required to explain the growth 
of the red sequence. The above estimates of growth 
through in situ star formation merely serve as a cau- 
tionary statement to indicate that a factor of 2 evolution 
in the red-galaxy mass density is rather modest and can 
be accounted for in various ways. 



Of course, if, for some reason, current measurements 
overestimate the evolution in the red-galaxy population, 
our results can be regarded as evidence for passive evo- 
lution of red galaxies, in terms of both star formation 
and mass assembly. We stress, however, that this is 
most likely not the case, given the mounting evidence 
for considerable evolution in the stellar mass density of 
red galaxies since z = 1. 

4.3. Comparison with Other Studies 

Numerous studies have analyzed the morphological 
and structural properties of galaxies in the local universe 
and at high redshifts. Even though it is not feasible to 
discuss all previous work, we put our work on the local 
galaxy population in the context of several other studies 
and compare our results with those from several recent 
efforts with the similar objective to measure the evolution 
with redshift of the early-type fraction outside massive 
cl usters. 

iGoto et all ()2003[ ) derive the MDR for luminosity- 
selected galaxies at 0.05 < z < 0.10. Their morphological 
classifications are simply based on the concentration of 
the light profile, and they find early-type galaxy fractions 
of ^ 40% — 60% at low densities, depending on the strict- 
ness of the classification crit erion for ea r ly typ es. This is 
consistent with our resuhs. iGoto et al.l ()2003l ) note that 
at densities below E ~ 1 Mpc~^ the early-type fraction 
changes less rapidly with density than at higher densities. 
This trend is even more pronounced in our mass-selected 
sample: at densities S < 1 Mpc^^ the early-type frac- 
tion is constant at ^ 40% and only increases at higher 
densities (see Fig. [10] and § 13. 2|) . 

The contributions of galaxies with different morpholo- 
gie s to the total stel lar mass density have been measured 
bv lBell et aP (|200^. who find that nowadays 76% of the 
stars in galaxies more massive than M — A x 10^" Mq 
reside in early-type galaxies. In our z 0.03 sample, 
we find a much lower value of 47%. The cause of this 
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difference is that Bell et al. use a single parameter (the 
concentration index) as a proxy for morphology, which 
correspon ds with visual classi fications 70% of the time 
(see also IStrateva et ahl |2001[ ). Our method uses the 
bumpiness parameter B in addition to the PSF-corrected 
concentration of the light profile, quantified by the Sersic 
parameter n. Thus, early-type galaxies are required to 
be not only highly concentrated but also smooth. This 
extension increases the agreement with visual classifica- 
tions to '^90% (see 12. 3p . This improvement is not due 
to smaller random errors: the contribution to the to- 
tal stellar mass density of highly concentrated galaxies 
in our z ~ 0.03 sample (those wi th n > 2.5 ) is 70 %, in 
good agreement with the result bv lBell et al.l (|2003D . and 
higher than that of early-type galaxies as defined by both 
visual classifications and by our B — n method. Many 
of the galaxies with n > 2.5 have late- type morpholo- 
gies, which is only revealed through visual classifications 
or the use of an additional p arameter such as bumpi- 
ness. As was demonstrated bv iKauffmann et al.l (|2004j ). 
the concentration distribution of the galaxy population 
provides crucial insight into the formation and evolution 
of galaxies. However, this structural property should be 
distinguished from morphology. The latter is related to 
the smoothness of the light profile in addition to its con- 
centration. It may well be that structure and morphology 
are correlated but distinct physical parameters that are 
affected by, e.g., galaxy mass and environment in differ- 
ent ways. 

Recently, several studies have addressed the evolution 
with redshift of the morphological mix of field galaxies. 
[Abraham et al.l (|2007f ) find that for galaxies more mas- 
sive than M = 4 X 10^° Mq the fraction of the stellar 
mass density residing in early-typ e galaxies is ~80 % at 
z 1, similar to the numbers from^ Bell et al.l (|2003f ) and 
our n-selected sample at low redshift. Abraham et al. 
use two parameters to distinguish early types from late 
types: in addition to the Gini coefficient (which, for early 
types and spiral galaxies, closely corresponds to the con- 
centration index) they use asymmetry as a morphological 
indicator. However, the addition of asymmetry docs not 
substantially improve the agreement with visual classi- 
fications for_galaxies with type Sb or later (see, e.g., 
IConselic3l2003l ). The method used bv [Abraham et all 
()2007f ) therefore effectively classifies galaxies by concen- 
tration. In our z ^ 0.8 sample, highly concentrated 
galaxies (with n > 2.5) contribute 72% to the total stel- 
lar mass density, similar to what Abraham et al. find 
(~80%). The stellar mass in early-type galaxies as de- 
fined by our B — n method or visual cla ssifications is 
lower (57%), consistent with the results bv lBundv et al.l 
(|2005), which are also based on visual classifications. We 
conclude t hat the difference betw een our numbers and 
those from [Abraham et al.l (|2007t l are due to the differ- 
ence between the classifying methods. Note that the rel- 
ative early-type galaxy mass density does not evolve be- 
tween z 1 and the present, independent of the applied 
morphological c lassifi cation method. 

ICapak et al.l ()2007D . who use the Gini coefficient to 
determine morphologies, find a significant decrease with 
redshift (in the range 0.4 < z < 1.2) of the early- type 
fraction in high-density regions (with E > 100 Mpc-2). 
At lower densities (with S < 50 Mpc"-^) they find a con- 
stant early- type fraction between z = 0.4 and 1.2 that 



is consistent with the early-type fraction that we find. 
However, when they compare their results with the local 
early-type fraction, measured with different classification 
methods, they do find significant evolutio n, also at the 
lowest densities. The sample constructed bv'Capa k et al.l 
fi2007i) is selected by luminosity, and is therefore not 
directly comparable with our sample. As was demon- 
strated above, selecting early-type galaxies by Gini co- 
efficient or concentration index alone overestimates the 
number of early-type galaxies as compared to visual clas- 
sifications or our B — n method. This, combined with the 
difference between luminosity- and mass-selected sam- 
ples, may conspire to yield an early-type fraction similar 
to that we find for visually classified galaxies in a mass- 
selected sample. 

We conclude that our results are either consistent with 
previous work or can be explained by the differences 
in approach. We note that the GOODS/ ACS imag- 
ing that we use to determine galaxy morphologies at 
z ^ 0.8 is much deeper than then imaging used by 
ICapak et al.l (I2007D and effect ively as deep as the data 
used bv [Abraham et"al] (|2007D in terms of S/N for red 
galaxies at z ~ 1. Moreover, our study is the first of its 
kind to quantify galaxy morphologies in the local uni- 
verse and at higher redshifts in an internally consistent 
manner. 

5. SUMMARY 

In this paper we examine the early-type galaxy 
(E-t-SO) fraction and its evolution for stellar-mass- 
selected, volume-limited samples of galaxies at redshifts 
< z < 1. At low redshift we estimate stellar masses 
and determine morphologies for 2003 galaxies at redshifts 
0.02 < z < 0.045 in the SDSS, complete down to a mass 
of Af = 4 X lO^*' Mq, 60% of the mass of a typical (L,) 
galaxy. In addition, we construct a similar sample of 
207 galaxies in the CDF-S at redshifts 0.6 < z < 1.0, 
complete down to the same mass limit. The stellar mass 
estimates are based on the relation between color and 
M/L, which has been proved to work without significant 
systematic effects with respect to kinematic mass mea- 
surements at low and high redshift. Morphologies are 
determined with an automated method, based on the 
Sersic parameter and the bumpiness of the residual. We 
compare this method with the traditional visual classifi- 
cations and find no systematic differences when we split 
the sample into two morphological classes (Sp-f Irr and 
E-t-SO). We estimate the local surface density with an 
nth nearest neighbor method. The density for the galax- 
ies in our samples is typically S ~ 1 Mpc~^ and ranges 
from '^0.01 to 100 Mpc~^. The sample characteristics 
are summarized in Table [TJ 

We find that for galaxies with masses higher than 
M = 4 X 10^° Mq the early- type fraction in the field 
and group environment probed by our samples has not 
changed significantly between z = 0.8 and the present. 
The early-type fraction is 43% ± 3% for the z ~ 0.03 
sample, and 48% ± 7% for the z ~ 0.8 sample. When 
we combine this with the unchanging early-type fraction 
of > 80% in dense (> 100 Mpc~^) cluster environments 
over the same redshift range and down to the same mass 
limit (Paper I), we find that the MDR has not evolved 
significantly between z ~ 0.8 and the present over at 
least 3 orders of magnitude in density. 
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In both our z ~ 0.03 sample and our z ~ 0.8 sample, 
65% — 70% of the galaxies (see also Tabled]) are located 
on the red sequence, in agreement with previous work. 
In addition, 65% — 70% of the galaxies in both samples 
have Sersic indices n > 2.5. Our samples follow the well- 
constrained trends that the stellar mass density of red 
galaxies has increased and the global star formation ac- 
tivity has decreased between z ~ 1 and the present. In 
fact, star formation in red galaxies could contribute sig- 
nificantly to the growth in the stellar mass density of 
red galaxies without strong evolution in the early-type 
galaxy fraction. 

Future work will include a more detailed analysis of the 
morphological mix, e.g., the separation of the late types 
into irregular and spiral galaxies and the separation of 
the early types into E and SO galaxies. Furthermore, 
studies in larger fields with high-resolution imaging will 
be used to measure the shape of the MDR at low den- 
sities in more detail. In particular, it is important to 
study intermediately dense regions, i.e., the infalling re- 



gions around clusters. At those densities the morpholog- 
ical mix of the galaxy population appears to undergo the 
most rapid change as a function of density. A crucial step 
forward in understanding the role of the environment and 
the role of the internal properties of galaxies in shaping 
their stellar populations and morphologies will be to in- 
clude the evolution of the density itself in analyses such 
as those carried out in this paper. Another challenge 
will be to establish at what redshift the MDR emerges, 
whether this coincides with the emergence of the red se- 
quence, and how this relates to the buildup of the stellar 
mass function. 
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